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ABSTRACT 

The  concepts  of  hennltlan  preserving  and  nonnegativity  preserving  subspaces 
of  complex  square  matrices  are  introduced.  Characterizations  of  such  subspaces 
are  obtained  and  applications  are  discussed. 
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1.  INTRODUCTION 


We  introduce  the  concepts  of  hermltian  preserving  and  nonnegativity  preserving 
subspaces  of  complex  square  matrices.  These  are  special  cases  of  the  concepts 
of  hermitian  preserving  and  positive  semidefinlte  preserving  linear  transformations 
given  in  dePlllis  [5]  and  later  investigated  by  Hlll[6],  Bose  and  Mitra  [4], 

Polulkis  and  Hill[12]  and  Barker,  Hill  and  Haertel  [1].  In  the  next  section,  we 
give  several  characterizations  of  hermltian  preserving  and  nonnegativity  preserving 
subspaces.  In  section  3,  we  use  these  results  to  investigate  the  existence  and 
characterization  of  hermitian  and  nonnegative  definite  solutions  to  linear  equa¬ 
tions.  Some  applications  are  discussed  in  section  4.  An  Important  application  in 
statistics,  orginlally  investigated  by  Pukelsheim  [13]  and  later  by  the  author  [9], 
is  also  given. 

The  following  notations  are  used  in  the  paper:  CnXm  and  Rn*m  respectively 

denote  the  vector  spaces  of  complex  and  real  matrices  of  order  nxm.  For  a  complex 
* 

matrix  A,  A  denotes  its  complex  conjugate  transpose  and  A  ,  its  transpose. 

The  inner  product  of  two  complex  matrices  A  and  B  of  the  same  order  is  assumed  to 

* 

be  the  trace  inner  product;  i.e,  <A,B>  -  tr  B  A.  The  range  of  a  linear  transforma¬ 
tion  T  is  denoted  as  R(T).  For  a  hermitian  matrix  A,  A+  and  A_  denote  the  positive 
and  negative  parts  of  A  respectively.  A+  is  obtained  from  the  spectral  decomposition 
of  A  by  deleting  the  negative  eigenvalues  and  the  corresponding  projectors.  A_ 

Is  then  defined  as  A+-  A.  If  a  square  matrix  A  is  hermitian  nonnegative  definite, 
we  denote  A^O. 

2.  HERMITIAN  AND  NONNEGATIVITY  PRESERVING  SUBSPACBS 

Let  A  be  a  subspace  of  CnXT>  and  let  P  denote  the  orthogonal  projector  onto  A 
(where  orthogonality  is  with  respect  to  the  trace  inner  product) . 


DEFINITION  (1)  A  is  said  to  be  hermitian  preserving  if  P(A)  is  hermitian  for  every 
hermitian  matrix  A. 

(2)  A  is  said  to  be  nonnegativity  preserving  if  P(A)^0  for  every  A^O. 

He  notice  that  nonnegativity  preserving  subspaces  are  those  which  preserve  the 
Loewner  ordering,  l.e.  if  A  and  B  are  hermitian  matrices,  then  A  is  nonnegativity 
preserving  iff  A-B>_0  Implies  P(A)  -  P(B)^0. 

THEOREM  1.  A  is  hermitian  preserving  iff  anyone  of  the  following  equivalent  conditions 
holds : 

(1)  A  £  A  whenever  Ac  A 

(ii)  P(A*)  ■  P(A)*  for  every  Ac  C0**1 . 

Proof:  (i)  Let  A* e  A,  for  every  Ac  A.  Let  B  be  a  hermitian  matrix  and  suppose 

P(B)*P(B)*.  Write  P(B)  -  P.  Then  ||B-P  ||2  -  ||B  -  ||  2  + 1|~~  -  P  ||2  + 

tr  (~_  _  P)*CB  _  Z±L_)  +  tr  (B  -  -  P).  Using  the  fact  that  B  is 

1  P+P*  *  1  P+P*  A 

hermitian,  the  last  two  terms  simplify  to  ytr  (B - ^ — )(P-P  )  +  jtr  (B j — )(P  -P), 

2  p | p*  2 

which  is  zero.  Thus  we  get  the  contradiction  1 1 B*-P  1 1  >  ||b - ^ — II  •  Con~ 

*  * 

versely,  let  A  be  hermitian  preserving  and  let  Ac  A.  Suppose  A  i  A.  Since  AfA 

*  * 

is  hermitian,  P  *  P(AtA-)  is  also  hermitian.  Then||A - 5 — 1|  ■  1 1 A— P  | [  + 

I! P  -  ^-||2  +  tr  (A-P)  (P  -  ^n-)  +  tr  CP  “  ^-)*(A-P).  The  last  two  terms  are 

A4A*  I 

zeros  since  A-Pc  A  and  P  -  —5 — c  A  (the  orthogonal  complement  of  A),  Thus  we 
*  2  * 

get  ||  A  -  ^2 A'  11^  >  1 1  A-P  || 2 ,  which  is  a  contradiction  since  is  the  hermitian 

matrix  closest  to  A;  see  Marshall  and  Olkln  ([8],  p.  264). 

(ii)  If  P(A*)  ■  P(A)*  for  every  Ac  Cn*n,  it  follows  that  P(A)  is  hermitian 

*  * 

whenever  A  is.  To  prove  the  converse,  write  A  ■  —j - i  ■  v  2 (i  ■  /-I  ). 

A+A*  i(A-A*)  .  .  *  A4A  .  .  i(A«-A  ) 

Since  ■-  and  ■  ' '■  are  hermitian  matrices  and  since  A  *  -5 —  +  i  2 
2  2 

it  follows  that  P( A  )  *  P( A)  when  A  is  hermitian  preserving.  Q 


9 


The  'if*  part  in  part  (ii)  of  Theorem  1  also  follows  from  the  lemma  in  Hill 


([6],  p.  260). 

COROLLARY  1.  A  hermitian  preserving  subspace  has  a  basis  consisting  of  hermitian 
matrices. 


THEOREM  2.  Let  A^ ,  A2 , .  • . ,  Ag  form  an  orthonormal  basis  for  A.  Let  denote  the 


(jk)**1  element  of  A^  (h«  1,2, . . .  ,s;  j  ,k  -  1,2,..., n).  Then  A  is  hermitian  preserving 
iff  anyone  of  the  following  equivalent  conditions  holds: 


(1)  The  matrices  J1<V|  +  V>\  *nd  V)  '  1  SJk)Ah  are 

for  j  ,k-  1,2, « • ,  ,n . 


,-h  .  -h 


-h 


-h 


hermitian 


(ii)  £  A^  A  Ah  is  hermitian, 
h*l 


Proof:  It  is  easy  to  see  that  the  set  of  hermitian  matrices  of  order  n  can  be  ob¬ 


tained  by  taking  real  linear  combinations  (.linear  combinations  with  real  coefficients) 


2  tti 

of  the  following  n  hermitian  matrices:  (a)  the  matrix  E^  having  j  diagonal 


th 


element  1  and  zeros  elsewhere  (j  -  l,2,...,n)  (b)  the  matrix  having  (jk)  and 


.th 


(kj)  elements  1  and  zeros  elsorhere  (j  <  k;  j  ,k-  1,2, . . .  ,n)  and  (.c)  the  matrix 


G  .  having  (jk)th  element  i,  (kj)th  element  -i  and  zeros  elsewhere  (j  <  k;  j,k«l,2,... 

J  O 


We 


have  P(B)  ■  £  (tr  A^B)^.  In  the  place  of  B  if  we  substitute  the  matrices  E^  ^ , 


F  .  and  G  .  defined  above,  we  get  part  (i)  of  the  theorem.  It  can  be  shown  that 

J  J  o  a 

-h  h  v  h  -h 


part  (i)  of  the  theorem  is  equivalent  to  )  a,  .a  ■  ][  (J  ,k,u,v- 1,2 . n) , 

s  ^  h-1  uv  h-1  w 

which  is  equivalent  to  [  A^  I  being  hermitian,  □ 

h"l 

Part  (ii)  of  the  above  theorem  also  follows  from  Proposition  1.2  in  dePillis 


[5]  or  Theorem  1  in  Hill  [6). 

THEOREM  3.  A  is  nonnegativity  preserving  iff  (i)  A  is  hermitian  preserving  and 
(11)  A  £  A  for  every  hermitian  matrix  Afi  A, 


Proof:  Suppose  A  Is  nonnegativity  preserving.  If  B  is  hermitian,  we  can  write 
B«B+-B_  and  P(B)  -  P(B+)  -  P(B_) .  Since  P(B+)^0  and  P(B_)^0,  P(B)  is  hermitian, 
showing  that  A  is  hermitian  preserving.  Let  Ac  A  be  hermitian.  Then  we  have 
A+-  A_-  A-  P(A)«  P(A+)  -  P(A_).  Hence  ||A+||2  +  ||aJ|2  -  HP(A+)  ||2  +  ||P(A_)  ||2  - 
2tr  P(A+)P(A_)<_  |P(A+)  ||2  +  ||P(A_)  J  *  (since  P(A+)>_0  and  P(A_)^0).  From  this,  it 
follows  that  A+-  P(A+)  or  A+c  A.  Conversely,  suppose  (i)  and  (ii)  in  the  theorem 
hold.  Let  B>_0.  Then  C  *  P(B)  is  hermitian  and  C  ^  0  if  C  is  not  nonnegative  defi¬ 
nite.  Now,  ||b-C|2  -  ||B-C+-C_ I2  -  1|b-C+|2+  [|c_ I  z-2tr  (B-C+)C_  >  ||B-cJ|2(since 
C+C_»0  and  tr  BC+>_0).  This  is  a  contradiction  since  C+  cannot  be  closer  to  B  than 
the  projection  C.  □ 

The  proof  of  the  "if"  part  is  similar  to  the  proof  of  Lemma  2  in  Pukelsheim 
[13]  and  the  proof  of  the  "only  if"  part  is  similar  to  the  proof  of  Lemma  1  in  Mathew 
[9]. 

COROLLARY  2.  A  nonnegativity  preserving  subspace  has  a  basis  consisting  of  hermitian 
nonnegative  definite  matrices. 

s 

THEOREM  4.  Let  A^,A2»...,Ag  be  an  orthonormal  basis  for  A.  If  ) 

h- 

A  is  nonnegativity  preserving.  This  condition  is  also  necessary  if  s-  2  or  if 
.  * 

there  exists  a  nonsingular  matrix  Q  such  that  Q  AQ  is  diagonal  for  every  hermitian 
matrix  Ac  A. 

Proof:  Let  A.  (j «  1,2, . . . ,s)  be  as  given  in  the  theorem  and  let  A>^0.  Then  P(  A)  - 
8  *  ^  8  *  * 

l  (tr  A^A)A^  and  P(A)>_0  for  every  A>_0  iff  £  (a  A^a)A.  >_0  for  every  vector  a  or 
h-1  ~  h-1 

equivalently 


ft  S  ft 

(a  ft  I)  (  l  Ab  •  A^  ( a  a  I)  >  0 


(I) 


for  every  vector  a.  This  proves  the  sufficiency  of  the  condition  as  asserted  in  the 

8  * 

theorem.  Conversely,  suppose  A  is  nonnegativity  preserving.  Then  £  A.  ft  A.  is 

h-1 

hermitian,  in  view  of  Theorem  3(i)  and  Theorem  2(ii).  Since  A  now  has  a  basis  con- 


5 


sis ting  of  hermitian  matrices  (Corollary  1) ,  if  there  exists  a  matrix  Q  as  specified 

* 

in  the  theorem,  then  Q  AQ  is  diagonal  for  every  Ac  A.  It  now  follows  that  (I)  holds 
8  * 

for  every  a  iff  £  A.  fl  A.  >_0.  When  s  -  2,  in  view  of  Corollary  2,  let  A..  j^O,  A-^0 
h-1  n  1 

be  a  basis  for  A.  Then,  using  Theorem  6.2.3  in  [15],  we  see  that  there  exists  a 

*  * 

nonsingular  matrix  Q  satisfying  Q  AjQ  and  Q  A2Q  are  diagonal.  Arguing  as  before, 
the  proof  can  be  completed.  □ 

s  * 

We  observe  that  the  condition  J  A^  1  A^>  0  is  similar  to  the  condition  given 

h-1 

in  Corollary  2.2  in  dePillis  [5].  This  condition  is  not  in  general  necessary  for  A 
to  be  nonnegativity  preserving.  For  example,  if  A  is  the  vector  space  of  2x2  complex 

fl  o' 


matrices,  then  clearly,  A  is  nonnegativity  preserving.  The  matrices  A^ 


b  o1 

,  A  -  — 

f°  l! 

,  .  1 . 

and  A.  -  — 

b  1 

.0  1 

3  /2 

1  Oj 

4  JS 

s  * 

-1  0. 

0  0 


form  an  orthonormal  basis  for  A  and 


do  not  satisfy  the  condition  \  A^  ft  A^^O. 

h-1 

REMARK  1.  The  existence  of  a  nonsingular  matrix  Q  satisfying  the  requirement  in 
Theorem  4  could  be  verified  as  follows:  When  A  is  nonnegativity  preserving,  let 
B.,B7,...,B  be  a  basis  for  A  satisfying  B.f^O  (j  -  1,2,. . . ,s) .  From  Theorem  3(ii) 

and  the  proof  of  Theorem  l(ii),  it  is  easy  to  see  that  such  a  basis  can  be  obtained. 

* 

Then  there  exists  a  nonsingular  matrix  Q  satisfying  Q  B  Q  is  diagonal  for  j-  l,2,...,s 

s  ^ 

iff  B.BqB^  -  B^BqB  for  j,k-  l,2,...,s.  Here  B-  -  \  B  and  BQ  denotes  a  generalized 

J  J  j»l  J 

inverse  of  Bq.  For  a  proof  of  the  above  assertion,  we  refer  to  Bhimasankaram 

([3],  Corollary  2)  or  Rao  and  Mitra  ([15],  Theorem  6,5.2). 

Suppose  A  is  a  subspace  of  real  square  matrices  of  order  n,  i.e.  Ac  R1L  We 

shall  call  A,  a  symmetry  preserving  subspace  if  P( A)  is  symmetric  whenever  A  is 

symmetric  (the  inner  product  under  consideration  is  <A,B>-  tr  B’ A  for  A,B  real 

matrices).  "Nonnegativity  preserving"  is  defined  similarly.  We  now  give  character- 


.n*n 


lzatlons  of  a  symmetry  preserving  or  nonnegativity  preserving  subspace  A  of  R 


n*n 


THEOREM  5.  A  subspace  A  of  R  is  symmetry  preserving  iff  anyone  of  the  following 
equivalent  conditions  holds: 

(i)  A'  c  A  whenever  Ac  A 

(ii)  P(A')  -  P(A) *  for  every  Ac  RnX“  . 

Proof:  The  proofs  of  (i)  and  the  'if'  part  of  (.11)  are  similar  to  the  correspond¬ 
ing  proofs  of  (i)  and  (li)  in  Theorem  1.  We  shall  now  prove  the  'only  if'  part 

of  (ii).  Suppose  A  is  symmetry  preserving  and  let  A., A,,..., A  form  an  orthonormal 

s  1  1  s 

basis  for  A.  Then,  for  any  Ac  RnX°,  P(A)  -  £  (trA/A)A..  In  view  of  part  (i) 

h-1  11  " 

of  the  theorem,  it  follows  that  the  matrices  A^,...,A'  also  form  an  orthonormal 

s  1  8 

basis  for  A.  Hence  P(A’)  -  j  (trA^A^A'^  «  P( A)'.  □ 


THEOREM 


Let  A  be  a  subspace  of  RnXn  and  let  A^.Aj^.-.A^  form  an  orthonormal 


basis  for  A.  Let  a^  denote  the  (jk)th  element  of  Afc  (h-  1,2,.  ...s;  j,k-  1,2,. . .  ,n) 
Then  A  is  symmetry  preserving  iff  anyone  of  the  following  equivalent  conditions 
holds: 

8  ^  ll  ll 

(i)  the  matrices  J  a  A.  and  J  (a,.  +  a  ,)A.  are  8y,m,,etrl-c  f°r  J»k“  l,2,...,n 

h-1  ^  h-1  •*  ^ 

8 

(ii)  £  (A^  fl  A^+  A^  a  A^l  is  symmetric, 
h— 1 

Proof :  The  matrices  E^  j-l,2,...,n  and  F  ^  j  <  k;  j  ,k-  1,2, , . .  ,n  defined  in  the 
proof  of  Theorem  2  form  a  basis  for  the  subspace  of  real  symmetric  matrices.  The 
proof  of  the  theorem  can  now  be  completed  along  the  lines  of  the  proof  of  Theorem  2 
THEOREM  7.  Let  A  be  a  subspace  of  RnXn.  Then 

(a)  A  is  nonnegativity  preserving  iff  (i)  A  is  symmetry  preserving  and  (ii)  A+c  A 
for  every  symmetric  matrix  Ac  A 

s 

(b)  If  A^,A2,...,Ag  are  as  in  Theorem  6  and  if  J  [ A^  a  A^+ A^  a  A^^O,  then  A  is 

“  -  h— 1 


nonnegativity  preserving. 


Proof:  The  rpoofs  of  (a)  and  (b)  are  similar  to  the  proofs  of  Theorem  3  and  the 
first  part  of  Theorem  4  respectively.  □ 

REMARK  A  is  a  subspace  of  Rn  n,  the  conclusions  in  Corollary  1  and  Corollary 

2  are  not  true  in  general  if  we  replace  "hermit ian"  by  "symmetric".  This  is  easily 
seen  by  taking  A  to  be  RnXn  itself. 

REMARKJJ.  When  A  is  a  subspace  of  real  symmetric  matrices,  all  the  conclusions  in 
Corollary  1,  Corollary  2,  Theorem  4  and  Remark  1  remain  valid  if  we  replace 
"hermitian”  by  "symmetric"  wherever  applicable. 


3.  HERMITIAN  AND  NONNEGATIVE  SOLVABILITY  OF  LINEAR  EQUATIONS. 

Let  T  be  a  linear  transformation  from  CnXn  into  C^Xt*.  In  this  section,  we  are 
interested  in  the  existence  and  characterization  of  hermitian  and  hermitian  nonnega¬ 
tive  definite  (nnd)  solutions  to  the  linear  equation  TX*  Q,  where  Qe  R(T)  (the 
range  of  T).  T*  denotes  the  adjoint  of  T  and  T+  denotes  the  maximal  generalized 
inverse  of  T  as  defined  in  Ben  Israel  and  Greville  ([2],  p.  318).  T+  is  also  the 
pseudoinverse  of  T  as  defined  in  Holmes  C17J,  pp,  216-226).  Then  R(T+)  ■  R(T  )  and 
T+T  is  the  orthogonal  projector  onto  RCT  )  (see  Theorem  2  on  p,  320  in  [ 2 J ) .  When 
T  is  a  matrix,  T+  coincides  with  the  Moore  Penrose  inverse  of  T.  We  say  that  TX*  Q 
is  hermitian  solvable  (nonnegat ively  solvable)  if  there  exists  a  hermitian  (respectivel 
hermitian  nnd)  matrix  X  satisfying  TX*  Q, 

DEFINITION ,  We  say  that  T+  characterizes  the  hermitian  solvability  (nonnegative 
solvability)  of  TX*  Q  if  the  existence  of  a  hermitian  (hermitian  nnd)  solution  to 

-f. 

TX*  Q  implies  that  T  Q  is  hermitian  (respectively  hermitian  nnd). 

Thus,  if  we  know  that  T+  characterizes  the  hermitian  solvability  (nonnega¬ 
tive  solvability)  of  TX*  Q,  in  order  to  examine  if  TX*  Q  has  a  hermitian  (hermitian 
nnd)  solution,  it  is  enough  to  check  if  T+Q  is  hermitian  (respectively  hermitian 


THEOREM  8.  (i)  For  every  Qe  R(T),  T+  characterizes  the  hermitlan  solvability  of 

TX*  Q  iff  R(T  )  is  hermitlan  preserving. 

(ii)  For  every  Qe  R(T) ,  T+  characterizes  the  nonnegative  solvability  of  TX- Q  iff 

* 

R(T  )  is  nonnegativity  preserving. 

*  + 

Proof ;  Suppose  R(T  )  is  hermitlan  preserving.  Then  T  TX  is  hermitlan  whenever  X 
is  hermitlan.  Thus  when  X  is  a  hermitlan  solution  to  TX-  Q,  T+TX-  T+Q  is  also  a 
hermitlan  solution  to  TX-  Q.  Conversely  let  T+  characterize  the  hermitlan  solva¬ 
bility  of  TX=  Q  for  every  Qe  R(T).  Let  Xq  be  any  hermitlan  matrix  and  let  TXq-  Qq. 
+  +  + 

Then  T  Qq  is  hermitlan.  But  T  Qq-  T  TXq  and  Xq  is  an  arbitrary  hermitlan  matrix. 

* 

Hence  R(T  )  is  hermitlan  preserving.  This  proves  (i).  (ii)  is  proved  similarly.  □ 
REMARK  4.  If  T  is  a  linear  transformation  from  RnXn  into  RpXq,  the  conclusions  in 
Theorem  8  remain  valid  if  we  replace  "hermitian"  by  "symmetric".  (Symmetric  solva¬ 
bility  of  TX*  Q  is  defined  analogous  to  hermitian  solvability.) 


4.  APPLICATIONS. 


p  xn  p  xp  s 

(a)  Let  A  e  C  J  and  B.c  C  J  J  for  j»l,2,...,s.  Let  \  p  *  p  and  consider 

J  j-1  J 

the  system  of  equations  AjXA^  "  ® j  (J  "  1,2 . s)  (assumed  to  be  consistent).  Here 

*  *  * 

TX  is  a  vector  consisting  of  the  s  matrices  (A.XA_,...,A  XA  )  and  R(T  )  consists  of 

3  *  PiXp1 

all  complex  wcrices  of  the  form  £  A.Z. Z  ,  where  Z .  e  C  J  .  From  Theorem  l(i), 

*  j-1  J  J  J 

if  follows  that  R(T  )  is  a  hermitian  preserving  subspace.  Hence  the  existence  of 
a  common  hermitian  solution  to  the  above  equations  can  be  verified  by  examining  if 
T+B  is  hermitian,  in  which  case  T+B  itself  gives  a  hermitian  solution  (B  denotes 
(B^ ,%2 »• •• »Bg)). 

(b)  Let  Aj  and  B^  be  as  defined  in  (a)  and  assume  further  that  Bj>_0  and  R(Bj)-R(Aj) 

*  n  *  *  * 

Let  A  be  a  matrix  satisfying  R(A  )  -  K(A^ :  A2:...:Ag).  Now  consider  the  equations 
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AjA  XAAj  ■  Bj,  j*  l,2,,..,s  and  suppose  we  want  to  examine  if  the  above  equations 
admit  a  common  hermitian  positive  definite  solution  X.  This  problem  arises  in 
connection  with  the  identification  problem  for  shorted  matrices,  see  the  proof  of 
Theorem  1  in  Mathew  and  Mitra  [10]  (see  also  Mitra  [11]  for  the  details  of  this 
problem)  .  Let  T  denote  the  corresponding  linear  transformation  (as  in  the  case  of 
(a)  above) .  We  notice  that  if  X  is  a  hermitian  nnd  solution  to  the  above  equations, 
then  P  XP  +  (I-P  )  is  a  hermitian  positive  definite  solution,  where  P  denotes  the 

A  A  A  A 

•k  .  A 

projection  matrix  A(A  A)  A  .  Thus,  it  is  enough  to  know  if  the  above  system  of 
equations  admits  a  common  hermitian  nnd  solution.  In  case  R(T  )  is  nonnegativity 
preserving,  this  can  be  verified  by  applying  Theorem  8(ii). 


(c)  Nonnegative  estimation  of  variance  components. 

Let  Y  be  a  random  Rn-vector  with  mean  vector  of  the  form  A0  and  dispersion 
k 

matrix  £  0  V  .  Here  A  is  a  known  matrix,  B  is  a  vector  of  unknown  parameters, 
j-1  J  J 

Vj  ( j  -  l,2,.,.,k)  are  known  real  symmetric  matrices  and  9^  ( j *  1,2,. . . ,k)  are  unknown 

parameters  known  as  variance  components.  0  »  (e^.O^. . . • ,0^)  is  such  that 
k  ^ 

£  0  V  ^0.  We  are  interested  in  estimating  the  linear  parametric  function 

j*lJ  J  + 

q’0  ■  q^0^  +  ’ ’ ,+^k0k’  **et  I-AA  and  let  T  be  the  linear  transformation  from 
Rnxn  into  Rk  defined  as  TX  -  (trXMV.M, , . . , trXMV  M) ' ,  for  Xe  Rn*n.  Here  Rn*n 

S  JL  K  S  S 

denotes  the  space  of  real  symmetric  matrices  of  order  n.  Then  it  is  known  that 


q’0  admits  a  translation  invariant  quadratic  unbiased  estimator  iff  TX  =  q  has  a 

solution  Xe  RgXn,  in  which  case  Y’MXMY  is  an  estimator  with  the  desired  properties 

(see  [14],  pp.  302-305).  When  qe  R(T),  q'0  admits  a  nonnegative  definite  quadratic 

unbiased  estimator  iff  TX»  q  admits  a  solution  Xf^O.  When  R(T  )  is  nonnegativity 

+  * 

preserving,  this  can  be  verified  by  examining  if  T  q>^0.  We  notice  that  R(T  )  is 
the  vector  space  spanned  by  the  matrices  MV^M, . . . ,MV^M.  Also,  Y’(T  q)Y  is  the  well 
known  minimum  norm  quadratic  unbiased  estimator  (MINQUE)  of  q'0,  introduced  by 


LO 


C.R.  Rao  ([14],  pp.  302-305).  Thus,  when  R(T  )  Is  nonnegativity  preserving,  the 
existence  of  an  nnd  quadratic  unbiased  estimator  of  q'0  is  equivalent  to  the  non¬ 
negativity  of  its  MINQUE.  Further  details  regarding  this  observation  are  given 
in  [13]  and  [9]. 

(d)  Nonnegativity  preserving  subspaces  of  Rn 

Let  S  be  a  subspace  of  Rn  and  let  P  denote  the  orthogonal  projector  onto  S. 
(Here  we  assume  that  orthogonality  is  with  respect  to  the  inner  product  (x,y)-y'x. 
More  general  inner  product  can  be  considered  and  the  modifications  in  such  a  case 
will  be  obvious) .  We  say  that  S  is  nonnegativity  preserving  if  Px  has  nonnegative 
components  whenever  the  vector  x  has  nonnegative  components.  Clearly,  S  is  nonnega- 

J^X|J 

tivity  preserving  iff  the  matrix  P  has  nonnegative  entries.  Let  Ae  R  .  Then, 
for  every  be  R(A)  ,  A+b  characterizes  the  nonnegative  solvability  of  Ax- b  iff  R(A') 
is  nonnegativity  preserving.  Thus,  when  R(A')  is  nonnegativity  preserving,  in  order 
to  verify  if  Ax-b,  be  R(A),  has  a  solution  x  with  nonnegative  components,  it  is 
enough  to  examine  if  the  vector  A  b  has  nonnegative  components,  in  which  case 
x*  A+b  is  a  desired  solution. 
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